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ABSTRACT 

We attempt to put constraints on different cosmological and biasing models by 
combining the recent clustering results of X-ray sources in the local (z < 0.1) 
and distant universe (z 1). To this end we compare the measured angular 
correlation function for bright (Akylas et al. 2000) and faint (Vikhlinin & 
Forman 1995) ROSAT^-iay sources respectively with those expected in three 
spatially flat cosmological models. Taking into account the different functional 
forms of the bias evolution, we find that there are two cosmological models 
which performs well the data. In particular, low-f2o cosmological models (J7a — 
1 — rio = 0.7) which contain either (i) high a™^^^ ~ 1.13 value with galaxy 
merging bias, b{z) oc (1 + z)^-^ or (ii) low a^^'^^^ — 0.9 with non-bias, h{z) = 1 
best reproduce the AGN clustering results. While rCDM models with different 
bias behaviour are ruled out at a high significance level. 

Keywords: galaxies: clustering- X-ray sources - cosmology: theory - large- 
scale structure of universe 



1 INTRODUCTION 

The study of the distribution of matter on large scales, 
based on diflterent extragalactic objects, provides impor- 
tant constraints on models of cosmic structure forma- 
tion. In particular Active Galactic Nuclei (AGN) can be 
detected up to very high redshifts and therefore provide 
information on how the X-ray selected sources trace the 
underlying mass distribution as well as the evolution of 
large scale structure (cf. Hartwick & Schade 1989). 

However, a serious problem here is how the lumi- 
nous matter traces the underlying mass distribution. 
Many authors have claimed that the large scale cluster- 
ing pattern of diflerent mass tracers (galaxies or clus- 
ters) is characterized by a bias picture (cf. Kaiser 1984). 
In this framework, biasing is assumed to be statistical 
in nature; galaxies and clusters are identified as high 
peaks of an underlying, initially Gaussian, random den- 
sity field. Biasing of galaxies with respect to the dark 
matter distribution was also found to be an essential in- 
gredient of CDM models of galaxy formation in order to 
reproduce the observed galaxy distribution (cf. Benson 
et al. 2000). Furthermore, diflerent studies have shown 
that the bias factor, b{z), is a monotonically increasing 
ftmction of redshift (cf. Fry 1996; Mo & White 1996; 
Matarrese et al. 1997; Moscardini et al. 1998; Tegmark 
& Peebles 1998; Basilakos & Plionis 2001). For example, 
Steidel et al. (1998) confirmed that the Lyman-break 
galaxies are very strongly biased tracers of mass and 
they found that b{z = 3.4)^ 6,4, 2, for SCDM, ACDM 
(fio = 0.3) and OCDM (Ho = 0.2), respectively. 



Studies based on the traditional indicators of clus- 
tering, like the two point correlation function, have 
been utilized in order to describe the AGN clustering 
properties. Our knowledge regarding the AGN cluster- 
ing comes mostly from optical surveys for QSO's (cf. 
Shanks & Boyle 1994; Groom & Shanks 1996; La Franca 
et al. 1998 and reference therein). It has been estab- 
lished, from Groom & Shanks (1996), that QSO's have 
a clustering length of ro = 5.4 ± l.l/i~^Mpc (with 
mean redshift of 1.27), while La Franca et al. (1998), 
analysing a sample of ~ 700 quasars in the redshift 
range 0.3 < z < 3.2, found ro = 6.2 ± 1.6^"^Mpc. 
It is very important to note that comparison of these 
clustering results in different redshifts rather favours a 
comoving model for the evolution of clustering. 

Similarly, Vikhlinin & Forman (1995) studied the 
angular clustering properties using a set of deep ROSAT 
observations. Carrera et al. (1998) combined two soft 
X-ray surveys (235 AGN), the ROSAT Deep Survey 
(Georgantopoulos et al. 1996) and the RIXOS survey 
(Mason et al. 2000) and found a spatial correlation 
length 1.5/i~^Mpc < ro < 5.5/i~^Mpc, depending on 
the adopted model of clustering evolution. Boyle & Mo 
(1993) analysing the EMSS survey which contains 183 
low redshift AGNs, found a feeble clustering signal on 
scales < lOft^^Mpc. Recently, Akylas et al (2000) us- 
ing 2096 sources detected in the ROSAT All Sky Sur- 
vey Bright Source Gatalogue (RASSBSG), derived the 
AGN angular correlation function in the nearby Uni- 
verse and utilizing Limber's equation obtained ro = 
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6.7 ± 1.0/i~^Mpc, assuming comoving clustering evolu- 
tion. 

In this paper we present the standard theoretical 
approach to estimate the angular correlation function 
w{9), using different models for the bias evolution in 
different spatially ffat cosmological models. Comparing 
the latter with observational results, we attempt to put 
constraints on the different cosmological and bias mod- 
els. The plan of this paper is the following: In section 2, 
(a) we present the calculation of the theoretical predic- 
tions for w{9), (b) we describe the models for bias evo- 
lution and (c) we discuss the AGN selection function. 
In section 3 we present the CDM spatially flat cosmolo- 
gies, while in section 4 we present the predictions for 
■w{9). The observational results are compared with the 
predictions in section 5 and finally in section 6 we draw 
our conclusions. 



2 THE INTEGRAL EQUATION 

For the purpose of this study we will utilize the rela- 
tion between the angular 'w{6) and spatial ^(r, z) two 
point correlation functions (cf. Maglioccheti et al. 1999 
and references theirin). As it is well known, this connec- 
tion can be done using the Limber equation (cf. Peebles 
1980). For example, in the case of a spatially flat Uni- 
verse (Slo-l-f2A = 1), the Limber equation can be written 
as 



w{e) = 2 



(1) 



where 0(a;) is the selection function (the probability 
that a source at a distance x is detected in the sur- 
vey) and X is the comoving coordinate related to the 
redshift through 



Ay 

Ho J, E{y) ' 
with 

E{z) = [no{i + zf + 



(2) 



(3) 



(see Peebles 1993). The mean surface density, M, on a 
survey of solid angle Q,s is: 



J\f = x^(j){x)dx = 



N{z)dz 



(4) 



where N(z) is the number of objects in the given sur- 
vey within the shell {z,z + dz). Therefore, combining 
the above system of equations, the expression for w{6) 
satisfies the form 

fr N^{z)E(z)dz f°° Ur,z)du 

^ ' ^ {jrN{z)dzr ^ ' 

The physical separation between two sources, separated 
by an angle 9 considering a small angle approximation 
is given by: 



(1 + 2) 



/ 2 , 2a2\2 
[U + X \ 



(6) 



Extending this picture, we quantify the evolution of 
clustering with epoch presenting the spatial correlation 
function of the X-ray sources as 



5(r, z) = Cmass(r)i?(z) 

with 

R{z) = D^{z)h\z) 



(7) 



(8) 



where D{z) is the linear growth rate of clustering (cf. 
Peebles 1993) Q being given by 



D{z)^ 



5noE{z) 



+ 

EHy) 



dy 



(9) 



and b{z) is the evolution of bias. On large scales, we 
utilize the standard procedure based on the linear the- 
ory predictions from the power spectrum P{k). On 
small scales we will follow the notation of Maglioc- 
cheti et al. (1999), who used the APM correlation func- 
tion (Maddox et al. 1990) in order to estimate the an- 
gular correlation function of radio sources. In partcu- 
lar, for F = 0.2 the above extrapolation is used for 
r < r"^* = 5.4/i~^Mpc while for other values of F we 
have re-scaled ^(r) using ^(r) cx 5(rFi/F2). 

2.1 Bias Evolution 

The concept of biasing between different classes of ex- 
tragalactic objects and the background matter distribu- 
tion was put forward by Kaiser (1984) and Bardeen et 
al. (1986) in order to explain the higher amplitude of 
the 2-point correlation function of clusters of galaxies 
with respect to that of galaxies themselves. 

The deterministic and linear nature of biasing has 
been chaUenged (cf. Bagla 1998; Dekel & Lahav 1999) 
and indeed on small scales (r < 10/i~^Mpc) there are 
significant deviations from b{r) = const. Despite the 
above, the linear biasing assumption is still a useful 
first order approximation which, due to its simplicity, 
it is used in most studies of large scale clustering (cf. 
Magliocchetti et al. 1999). In this paper however, we 
will not indulge any variation of bias with scale but 
rather, working within the paradigm of linear and scale- 
independent bias. Therefore, we shortly describe some 
of the bias evolution models in order to introduce them 
to our analysis (see eq. |^ and eq. 

• No evolution of bias (BO): This model considers 
constant bias at all epochs: 



b{z) = bo ~ bT 



To 



7/2 



(10) 



where r"^' — 5.4/i~^Mpc is the correlation length in co- 
moving coordinates estimated by the APM correlation 
function (Maddox et al. 1990) and ro is the correspond- 
ing length for the X-ray sources. Finally, b^^^ — 1/(7™^°° 
is the present bias of optical galaxies relative to the dis- 
tribution of mass and ag^'^^ is the mass rms fluctuations 
in sphere of radius 8/i~^Mpc. Using the above ideas, if 
ones assumes that ro > r"^' , then it is quite obvious that 

* D{z) = (1 -I- z)~^ for an Einstein-de Sitter Universe. 
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the X-ray sources are indeed more biased with respect 
to optical galaxies by the factor (ro/r"''*)^''^. In case 
that h{z) = 1, we have the so called non-bias modelQ. 

• Test Particle or Galaxy Conserving Bias (Bl): This 
model, proposed by Nusser & Davis (1994), Fry (1996), 
Tegmark & Peebles (1998), predicts the evolution of 
bias, independent of the mass and the origin of halos, 
assuming only that the test particles fluctuation field is 
related proportionally to that of the underlying mass. 
Thus, the bias factor as a function of redshift can be 
written: 



h(z) = 1 + 



jbo - 1) 
D{z) 



(11) 



where 6o is the bias factor at the present time. Bagla 
(1998) has found, that in the range < z < 1, the 
above formula describes well the evolution of bias. 

• Merging Bias Model (B2): Mo & White (1996) have 
developed a model for the evolution of the correlation 
bias, using the Press-Schechter formalism. Utilizing a 
similar formalism, Matarrese et al. (1997) extended the 
Mo & White results to include the effects of different 
mass scales (see also Moscardini et al. 1997; Bagla 1998; 
Catelan et al. 1998; Magliocchetti et al. 2000). In this 
case we have that 



(12) 



with P ~ 1.8. 



2.2 Selection Function 

In flux-limited samples, it is well known that there is a 
degradation of sampling as a function of distance from 
the observer (codified by the so called selection func- 
tion). The latter also depends on the evolution of the 
source luminosity function. Thus for the X-ray sources 
the selection function can be written as 



$(L,,z)dL . 



(13) 



In this work we used a luminosity function of the form 
assumed by Boyle et al. (1993), which takes into account 
the cosmological evolution of the QSO's in the form of 
pure luminosity evolution. 



3 COLD DARK MATTER (CDM) 
COSMOLOGIES 

In this section, we present the cosmological models that 
we use in this work. For the power spectrum of our 
CDM models, we consider P(k) ^ k"T'^{k) with scale- 
invariant (n — 1) primeval inflationary fluctuations. 
We utilize the transfer function parameterization as in 
Bardeen et al. (1986), with the corrections given ap- 
proximately by Sugiyama's (1995) formula: 



T{k) 



ln(l + 2.34g) 
2.34q 



l + 3.89q + {W.lqy + 



{5A6qf + (6.71g)*] 



-1/4 



with 



k 



(14) 



where = 27r/A is the wavenumber in units of h Mpc~^ 
and ili, is the baryon density. In this analysis, we have 
taken into account three different cold dark matter mod- 
els (CDM) in order to isolate the effects of different pa- 
rameters on the X-ray sources clustering predictions. 

The rCDM0 and ACDM (see Martini & Wein- 
berg 2000) are approximately COBE normalized and 
the latter cosmological model is consistent with the 
results from Type la supernovae (Riess et al. 1998; 
Perlmutter et al. 1999). In the same framework, the 
rCDM and ACDM models have F ~ 0.2, in approx- 
imate agreement with the shape parameter estimated 
from galaxy surveys (cf. Maddox et al. 1990; Peacock & 
Dodds 1994) and they have fluctuation amplitude in 8 
ft^^Mpc scale, (t™^°°, consistent with the cluster abun- 
dance, = 0.55^-° '^ (Eke, Cole, & Frenk 1996). 

Furthermore, in order to investigate cosmological 
models with high value of ag^'^"^ , we include a new model 
named ACDM2 (cf. Cole et al. 1998). The crg^^^" value 
for the latter cosmological model is in good agreement 
with both cluster and the 4-years COBE data with a 
shape parameter F — 0.25. Therefore, it is quite obvi- 
ous that two of our models (rCDM and ACDM) have 
the same power spectrum and geometry but different 
values of Qo and ag^^^^, while the two spatially flat, 
low-density CDM models (ACDM and ACDM2) have 
different ag and F respectively. In table 1, we present 
the normalizations for the above speciflc cosmological 
models as well as the values of the AGN bias, bo, at the 
present time (described by eq. Thus it turns out 
that in ACDM and rCDM models the distribution of 
X-ray sources is "biased" relative to the distribution of 
mass; while the ACDM2 model is almost non "biased" . 



4 THEORETICAL PREDICTIONS 

In this section we will flrst present results on the pre- 
dicted angular correlation function w{9) (see Figure 1) 
estimated for three different CDM spatially flat cosmo- 
logical models taking into account the functional forms 
of the bias evolution (BO, Bl, B2) introduced in the 
previous section. From figure 1, it is obvious that the 
amplitude of w{0) increases for models with bias evolu- 
tion in all cosmological models. This is to be expected 
due to the fact that the amplitude of the angular corre- 
lation function is affected by the expression R{z), espe- 
cially at high redshifts. Indeed, we observe that models 



t From now on we consider b{z) = 1 as a (BO) bias model. 



t This model could correspond to the decaying neutrino 
model. 



4 Basilakos, S. 



Table 1. Small Scale Normalizations 



Model 


Ho 


Ha 


h 


mass 


bo 




rCDM 


1.0 


0.0 


0.50 


0.55 


2.18 


0.050 


ACDM 


0.3 


0.7 


0.65 


0.90 


1.33 


0.036 


ACDM2 


0.3 


0.7 


0.60 


1.13 


1.06 


0.035 




10 



^ 1 
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10 



10 
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Figure 1. Theoretical predictions of the angular correlation 
function w{9). 



0.1 1 
z 

Figure 2. The function R{z) = D^{z)b'^(z) as a function of 
redshift, for different bias evolution models. 



with increasing bias as a function of redshift, give a 
stronger clustering signal (in small angular scales), rel- 
ative to models with non-bias (BO). For the latter bias 
behaviour, it is interesting to see, that the predictions 
for w{9) for the case of the rCDM and ACDM2 are 
almost the same. 

Therefore, in order to understand better the effects 
of AGN clustering, we present in figure 2 the quantity 
R{z) — D^{z)h^{z) as a function of redshift for the three 
cosmological models, utilizing at the same time different 
bias evolution. It is quite obvious that the behaviour of 
the function R{z) characterizes the clustering evolution 
with epoch. Figure 2, for example, clearly shows that the 
bias at high redshifts has different values in different 
cosmological models. In particular for the high a"^^^^ 
low-f2o flat model (ACDM2) the distribution of X-ray 
sources is only weakly biased, as opposed to the strongly 
biased distribution in the rCDM cosmological model. 

Indeed the different functional forms of b{z), pro- 
vide clustering models where: 

• AGN clustering is a decreasing function with red- 
shift for (BO), 

• AGN clustering is roughly constant for (Bl). How- 
ever, the ACDM2-B1 model gives lower R{z) simply be- 
cause the higher a"^^^^ normalization largely removes 
the clustering difference between the two other flat 
cosmological models with low cr™'"'° normalizations. In 
other words, the present bias value of the above model 



is almost ~ 1, which gives clustering behaviour similar 
to the ACDM2-B0. 

• AGN clustering is a monotonically increasing func- 
tion of redshift for (B2). 



5 APPLICATION TO THE DATA 

In Figure 3 we compare the angular correlation func- 
tion of a sample of 2096 sources with a total sky cov- 
erage of 4.9sr detected in the ROSAT All-Sky Sur- 
vey Bright Source Catalogue (see Akylas et al. 2000) 
with that predicted in various flat cosmological models. 
Considering a two point angular correlation function of 
the form w{9) = {9 /9oY~'' , the above authors found 
9o = 0.062°, 7 = 1.8 and spatial correlation length of 
ro ~ 6.5 ± l.Oft-^Mpc and To ~ 6.7 ± 1.0/i~^Mpc for 
stable and comoving clustering evolution respectively, 
similar to the optically selected AGN. Due to the fact 
that the above estimates have been focused on the lo- 
cal X-ray Universe, in this work we use complementary 
observational results from Vikhlinin & Forman (1995). 
They have analysed a set of deep 7? 05^ T observations 
with a total sky coverage of 40 deg^, in order to inves- 
tigate the clustering properties of faint X-ray sources. 
Therefore, taking into account the correction for the 
amplification bias they claimed that the two point an- 
gular correlation function is well described by a power 
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Figure 3. Comparison of the predicted angular correlation 
function for various cosmological models witli tliat of the lo- 
cal AGN distribution, estimated by Akylas et al. (2000). Er- 
rorbars are determined by assuming Poisson statistics. The 
continuous dot-dash line represent the best fit to Wx{d) de- 
rived by the above authors. 



Table 2. probabilities ("P^ 2) °f consistency between 
RASSBSC data and models. 



Comparison Pair 


■pB 


RASSBSC 


- ACDM-BO 


0.81 


RASSBSC 


- ACDM-Bl 


0.067 


RASSBSC 


- ACDM-B2 


0.019 


RASSBSC - 


ACDM2-B0 


0.19 


RASSBSC - 


ACDM2-B1 


0.42 


RASSBSC - 


ACDM2-B2 


0.57 


RASSBSC 


- rCDM-BO 


0.19 


RASSBSC 


- rCDM-Bl 


2.88x10 


RASSBSC 


- rCDM-B2 


3.00x10 



low with 7 = 1.8. Thus, in figure 4 we plot their results 
and the estimated angular correlation function for all 
nine models. 

In order to quantify the differences between mod- 
els and data, we perform a standard test, for the 
bright (RASSBSC) and faint sources respectively, and 

we present the ^^^2!] results in tables 2 and 3. Com- 
paring the statistical results for both (a) high (faint 
sources) and (b) low (bright) redshift regimes, we can 
point out that for the bright X-ray sources the only 
models that are excluded by the data, at a relatively 
high significance level, are ACDM-B2, tCDM-B1 and 



8 Where N = B,F for bright (RASSBSC) and faint X-ray 
sources respectively. 




Figure 4. Comparison of the predicted angular correlation 
function for various models with that estimated by Vikhlinin 
& Forman (1995). The continuous dot-dash line represent the 
best fit to Wx (0) derived by the above authors. 



Table 3. probabilities {'P^^2 ) °f consistency between 
faint ROSAT data and models. 



Comparison Pair 


■pF 




faint ROSAT 


- ACDM-BO 


0.037 




faint ROSAT 


- ACDM-Bl 


3.67 xlO 


-4 


faint ROSAT 


- ACDM-B2 


1.72 xlO 


-6 


faint ROSAT - 


ACDM2-B0 


3.44 XlO 


-4 


faint ROSAT - 


ACDM2-B1 


3.21 XlO 


-3 


faint ROSAT - 


ACDM2-B2 


0.036 




faint ROSAT 


- tCDM-BO 


4.64 XlO 


-3 


faint ROSAT 


- rCDM-Bl 


1.07x10" 


12 


faint ROSAT 


- TCDM-B2 


4.73 xlO" 


-15 



rCDM-B2. Interestingly, for the faint objects the ex- 
cluded models are ACDM-Bl, ACDM-B2, ACDM2-B0, 
ACDM2-B1, rCDM-BO, rCDM-Bl and rCDM-B2. The 
above differences between the two kind of populations 
are to be expected simply because the cosmological evo- 
lution plays an important role on large scale structure 
clustering due to the fact that the high redshift objects 
are more biased tracers of the underlying matter dis- 
tribution with respect to the low redshift objects (cf. 
Steidel et al. 1998). Also, from the faint X-ray sources 
results we would like to point out that there is not a 
single rCDM model that fits the data. 

If we make the reasonable assumption that there 
is no correlation between the two X-ray populations, 
mostly due to the large distances involved, we can con- 
sider the RASSBSC and faint ROSAT catalogues as 
being independent of each other. Under this assump- 
tion the previous statistical tests can also be consid- 
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Table 4. The joint probabilities P. 



Fitted Models 


P 




ACDM-BO 


0.14 




ACDM-Bl 


2.86 xlO" 


-4 


ACDM-B2 


5.96 xlO 


-7 


ACDM2-B0 


7.00 xlO 


-4 


ACDM2-B1 


0.01 




ACDM2-B2 


0.10 




tCDM-BO 


7.08 XlO" 


-3 


rCDM-Bl 


1.41x10- 


18 


rCDM-B2 


7.29x10" 


22 



ered as independent. Indeed performing a standard 
Kolmogorov-Smirnov test comparing the two different 
spatial correlation functions, described by Akylas et al. 
(2000) and Vikhlinin & Forman (1995) wc find that they 
are significantly different. In this framework the joint 
(overall) probability can be given by the following ex- 
pression: 

P = V^^.Vl^.[l-HV^,^.Vl^.)\ . (15) 

In table 4 we can see the corresponding joint probabili- 
ties for all our models. This overall statistical test proves 
that the ACDM-BO and ACDM2-B2 models fit well the 
observational data at a relatively high significance level. 

Note that we have tested the AGN clustering pre- 
dictions using also the SCDM (COBE normalized from 
Bunn & White 1997, with h = 0.5, V = 0.5 and 
^mass _ -^ 22) cosmological model and wc found that 
it is excluded by the data at a high significance level 
^y,'^ — 10^*' — 10^^). Finally, a possible contamina- 
tion of the AGN samples by stars, which can be as- 
sumed at first order to be distributed randomly on the 
sky, will lower the true amplitude of the AGN wiO). Wc 
have crudely tested the effect on our model comparion 
results of an AGN amplitude drop of ~ 10% and found 
no qualitative differences whatsoever. 

We should conclude that the behaviour of the ob- 
served angular correlation function of the X-ray sources 
is sensitive to the different cosmologies but at the same 
time there is a strong dependence on the bias models 
that we have considered in our analysis. By separating 
between low and high redshift regimes, we obtain re- 
sults being consistent with the hierarchical clustering 
scenario, in which the AGN's are strongly biased at all 
cosmic epochs (cf. Magliocchetti et al. 1999). 



6 CONCLUSIONS 

We have studied the clustering properties of the X-ray 

sources using the predicted angular correlation function 
for several cosmological models. We parametrize the 
predictions for 111(6) taking into account the behaviour 
of b{z) for a non-bias model, a galaxy conserving bias 
model with b[z) (x (1 + z) and for a galaxy merging 
bias model with h(z) cx (1 + z)^ *. Utilising the mea- 
sured angular correlation function, for faint and bright 
X-ray ROSAT sources, estimated in Vikhlinin & For- 
man (1995) and Akylas et al. (2000) respectively, we 



have compared them with the corresponding ones pre- 
dicted in three cosmological models, namely the rCDM, 
ACDM and ACDM2 (with a high c7^'^== value). We find 
that the models that best reproduce the observational 
results are: 

• ACDM2 model (f^A = 1 - f2o = 0.7) with high 
^mass _ -[^ g^j^j hi&s evolution described by b[z) oc 

• ACDM model (Qa = 1 - = 0.7) with low 
^mass ^ Q g ^j-gg evolution by b{z) = 1. 
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